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Abstract 



It is well known that solutions to the Fourier- 
Galerkin truncation of the inviscid Burgers equa- 
tion (and other hyperbolic conservation laws) do 
not converge to the physically relevant entropy so- 
lution after the formation of the first shock. This 
loss of convergence was recently studied in detail in 
[S. S. Ray et al, Phys. Rev. E 84, 016301 (2011)], 
and traced back to the appearance of a spatially lo- 
calized resonance phenomenon perturbing the solu- 
tion. In this work, we propose a way to remove this 
resonance by filtering a wavelet representation of 
the Galerkin-truncated equations. A method pre- 
viously developed with a complex-valued wavelet 
frame is applied and expanded to embrace the use 
of real- valued orthogonal wavelet basis, which we 
show to yield satisfactory results only under the 
condition of adding a safety zone in wavelet space. 
We also apply the complex-valued wavelet based 
method to the 2D Euler equation problem, show- 
ing that it is able to filter the resonances in this 
case as well. 



1 Introduction 

Due to the intrinsic limitations of computers, solv- 
ing a nonlinear partial differential equation numer- 
ically actually means solving its truncation to a fi- 
nite number of modes, where, in favorable cases, 
the truncated system closely approaches its con- 
tinuous counterpart. But sometimes the trunca- 
tion has drastic effects which completely destroy 
the desired approximation. The first historical ex- 
ample for which this happened was probably the 
symmetric finite difference scheme designed by von 
Neumann in the 1940s for nonlinear conservation 
laws. As recalled in [6 , it was indeed shown in 
the 1980s that, when applying this scheme even to 
the simplest case of the ID inviscid Burgers equa- 
tions, convergence to the correct solution is lost at 
the appearance of the first shock. Other schemes, 
specifically designed to dissipate kinetic energy at 
the location of shocks, do not suffer from this limi- 
tation and yield the desired solution. 

This matter of convergence was investigated in 
[T6] for another important scheme, namely Fourier- 
Galerkin truncation, where only the equations for 
Fourier modes with wavenumbers below a certain 
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cut-off are solved, the other modes being set to zero. 
Using the conservative character of the truncation 
and the nonlinear structure of the equations, the 
author was able to prove that even weak conver- 
gence to the physical solutions was ruled out once 
the latter started to be dissipative. This loss of 
convergence was scrutinized more closely in the re- 
cent work [12] , which showed that in the truncated 
system shocks become sources of waves that per- 
turb the numerical solution throughout its spatial 
domain. This is possible because Fourier- Galerkin 
truncation is a non-local operator in physical space, 
instantaneously removing all modes above the trun- 
cation wavenumber. Furthermore, these waves res- 
onantly interact with the flow at locations where 
the velocity is the same as their phase velocity, giv- 
ing rise to strong perturbations localized around 
these positions which eventually spread and cor- 
rupt the numerical solution. 

The aim of the present work is to show how the 
resonances can be eliminated by filtering the solu- 
tion in a wavelet basis, a possibility which was al- 
ready pointed out in [12] . The Burgers equation has 
been chosen as a toy model because its entropy solu- 
tions can be computed analytically, enabling direct 
comparison with numerical results. An important 
point to keep in mind though is that the analytical 
solutions are dissipative even in the inviscid limit, 
a phenomenon known as dissipative anomaly, while 
the Galerkin-truncated ones never dissipate energy 
if the viscosity is set to zero. Therefore a numerical 
solution can approach the exact solution only if it 
finds a way to dissipate energy, as is achieved by 
our method through the filtering process described 
further down. In fact, as discussed in [10] [9] and 
references therein, many filtering mechanisms are 
known empirically to achieve this task (see also the 
recent review in [5]). However, the precise effect of 
these filtering methods on the resonances shown by 
[T2] has not been fully clarified yet. 

To get insight into the formation of the resonance 
we start by performing a continuous wavelet analy- 
sis of the Galerkin-truncated solutions to the invis- 
cid Burgers equation. Such a representation unfolds 
the solution in both space and scale in a continu- 
ous fashion. It thus allows to visualize at which 
wavenumbers and positions the resonances are gen- 
erated and subsequently propagated. 

Afterwards, the wavelet filtering method analo- 
gous to Coherent Vorticity Simulation (CVS), al- 
ready proposed to solve Burgers equation [10] [9] , is 
applied here with the same initial conditions used 
in [12] . To demonstrate that the method is well 
suited for regularizing the solution, the equation 
is solved in Fourier space using a pseudo-spectral 
approach, but after each time step the solution is 



expanded over a frame of complex- valued wavelets, 
filtered with an iterative procedure introduced in 
PQ, and then reprojected onto the Fourier basis for 
computing the next time step. 

We then go further and propose the use of real- 
valued orthogonal wavelets instead of the redun- 
dant complex-valued wavelets. Since the former do 
not enjoy the translational invariance property of 
the latter, satisfactory solutions can only be ob- 
tained by keeping the neighbors of the retained co- 
efficients, i.e., adding a safety zone in wavelet coef- 
ficient space to account for the shocks translation 
and the small scale generation, a procedure success- 
fully applied in previous works for 2D and 3D flows 
[4] [15] [11] . The quality of the approximations ob- 
tained for the different filtering methods is assessed 
by computing a global error estimate. 

Finally, since [12 also discusses the presence of 
resonances in the Galerkin-truncated 2D incom- 
pressible Euler equations, we accordingly study the 
effect of the complex- valued wavelet method in this 
case. First results in that same direction can be 
found in [9]. 

2 ID inviscid Burgers equa- 
tion 

2.1 Continuous wavelet analysis 

Our starting point is the inviscid Burgers equation, 
written in conservative form 

d t u+^d x u 2 = 0, (1) 

u being velocity, t time and x space, plus peri- 
odic boundary conditions, and taking the same har- 
monic initial condition as in [12 (the domain size 
being normalized to 1): 

^o(^) = sin(27nr) + sin(47rx + 0.9) + sin(67rx). (2) 

In [12] the authors observed that, when solving the 
Galerkin-truncated version of ([!]) with a pseudo- 
spectral code, fine scale oscillations appear all over 
the solution right after the formation of the first 
singularity in the exact solution, followed by the 
emergence of two bulges around the points having 
the shock velocity with positive velocity gradient. 
These bulges then grow and start to perturb the 
solution, initiating the equipartition process pre- 
dicted by T.D. Lee [8]. As explained in [12 , the 
bulges are due to a resonant interaction between a 
truncation wave, excited by the Gibbs oscillations 
coming from the Galerkin truncation, and the lo- 
cations where the velocities are close to the phase 
velocity of the wave. 
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Figure 1: Base 10 logarithm of moduli of CWT coefficients of the Galerkin-truncated solution, at time 
t = (a), t = 0.02749 (b), t = 0.03505 (c), t = 0.03538 (d), t = 0.03648 (e), t = 0.03998 (f), 
t = 0.05897 (g), and t = 0.19989 (h). Each subfigure shows the solution (on the top) and below the base 
10 logarithm of the modulus of the corresponding CWT coefficients. The corresponding wavenumber 
spectrum is plotted vertically on the left. 



To follow the formation of resonances and the 
subsequent spreading of the fluctuations, let us first 
consider the continuous wavelet transform (CWT) 
of the numerical solution at different time instants. 
All computations were performed using a 4 th or- 
der Runge-Kutta time evolution scheme with St = 
0.1257V -1 as time step, up to a Galerkin truncation 
wavenumber K = 8192. For efficiency, the nonlin- 
ear term is computed pseudo-spectrally on a col- 
location grid having 3K points, which ensures full 
dealiasing. The CWT coefficients are calculated as 
the inner products of the velocity u{x) at a given 
instant t with a set of wavelet functions ^i yX >(x) of 
scales t centered around positions x' where for the 
mother wavelet we have chosen the complex- valued 
Morlet wavelet for its good analysis properties [3]. 
The results, presented in Fig.[l] show the logarithm 
of the modulus of wavelet coefficients at different 
positions x' and scales t (represented by the equiv- 
alent wavenumbers k = ^r, k^ being the centroid 
wavenumber of the chosen wavelet [33]). The hori- 
zontal black line indicates the Galerkin truncation 
frequency and the velocity fields themselves are also 
shown at the top of each figure for convenience. 



Figures [Ta| and [Tb] show respectively the har- 
monic initial condition and how the precursors of 
the shocks develop. Figure [Tc] shows the solution 
when the first preshock reaches the cut-off scale and 
becomes a shock, i.e., when non negligible energy 
reaches the scale indicated by the horizontal black 
line. We observe that the first resonances appear 
immediately after that (note the small time interval 
between Figs. [Tq and Id) and then spread all over 



space. Figure [Tel shows the formation of the bulges 
around the resonant locations. They stretch un- 
til they reach the Galerkin scale and then generate 
more truncation waves, as shown in Fig. [If] Af- 
ter that, perturbations at all scales start to spread 
throughout the solution, and even more so when 
the second shock is formed, as in Fig.[lg[ For much 
longer time the solution then becomes very noisy 
(Fig. lh), on its way towards equipartitiorj^] 



2.2 Elimination of resonances using 
complex-valued Kingslets 

As explained in [T2] . and as we have seen from the 
wavelet analysis of the previous section, the fail- 
ure of the Fourier- Galerkin scheme to reproduce the 
correct solution can be traced back to the amplifi- 
cation of truncation waves by a resonance mecha- 



1 Videos with the time evolution of the coefficients were 
made available on-line for the interested reader as sup- 
plementary material to this paper, and also at http: 
//www . youtube . com/watch?v=WX2YIHGR7LA]and |http : //www . | 
youtube . com/watch? v=j 4Vf BGgSy30 



nism. To suppress these resonances, a dissipation 
mechanism has to be introduced in the numerical 
scheme, in a way which does not affect the nonlin- 
ear dynamics. This procedure is sometimes called 
regularization of the solution. In this section, we 
show by numerical experiments how the resonances 
are suppressed by the CVS-filtering method, which 
was first applied to the inviscid Burgers equation 
in [10], and recall its interpretation in terms of de- 
noising. 

The algorithm proposed by [10 is as follows. 
Starting from the Fourier coefficients of the velocity 
field Uk for \k\ < K at t = t n : 

1. Time integration. The Fourier coefficients of 
the velocity field are advanced in time to t = 
t n+ i using the 4 th order Runge-Kutta scheme 
described in Sec. 12.11 



2. Inverse Fourier transform. The velocity field 
at t = t n+ \ is reconstructed from its Fourier 
coefficients on a grid with N = 2K points. 

3. Forward wavelet transform. The velocity field 
is written in wavelet space as 

J-l 2 j 

U(x) = (0| u) <t>{x) + ^ Yl ^3i\ U ) i>ji( x )> 
j=0 i=l 

(3) 

where ipji are the wavelet functions, <p the 
associated scaling function and the indexes 
j and i denote scale and position respec- 
tively. Each inner product, defined as (f\g) = 
Jo f( x )*9( x )dx, corresponds to a wavelet coef- 
ficient. 

4. Application of the CVS filter. The coefficients 
whose modulus are below a threshold T, so- 
called incoherent coefficients, are discarded, 
and T is determined at each time step in an 
iterative way following pQ. It is initialized as 
To = qy/E/N, q being a compression parame- 
ter and E being the total energy, then succes- 
sive filterings are made as T is recalculated in 
sub-step n + 1 as 



+i 



qo 



(4) 



until T n+ i = T n . Here are the wavelet co- 
efficients below the threshold T n and a[-] rep- 
resents the standard deviation of the set of co- 
efficients between brackets. 

Inverse wavelet transform. The coefficients 
above the final threshold represent the coher- 
ent part of the signal and are used as input to 
an inverse fast wavelet transform. 
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6. Forward Fourier transform. The Fourier coef- 
ficients of the filtered velocity field are com- 
puted, and the cycle can proceed onward. 

There are two choices left to be made in this 
algorithm: the wavelet basis used in steps 3 and 
5, and the parameter q in step 4. As shown in 
[9], this version of the algorithm performs badly if 
real- valued orthogonal wavelets are used, but works 
very well when using translation invariant complex- 
valued wavelets called Kingslets, introduced in [7] 
and first proposed in [10] for this application. Note 
that Kingslets were constructed to have almost van- 
ishing energy in the negative wavenumber range, 
which (as explained in [7j) implies that filtering in 
wavelet space is almost a translation invariant op- 
erator (i.e., it commutes with spatial translations 
of the signal). This is a desired feature for Burg- 
ers equation since shocks translate and cannot be 
properly tracked with a real-valued wavelet basis, 
whose coefficients are not stable enough due to the 
loss of translational invariance, giving poor filter- 
ing results. Therefore, we stick to this choice in 
this section, but we will discuss below how the al- 
gorithm can be modified to authorize other choices. 

Concerning the dimensionless number q in step 4 
of the algorithm, it controls the severity of the filter, 
since increasing q enlarges the set of discarded coef- 
ficients. Its value defines a certain balance between 
regularization and approximation quality, and also 
influences the compression rate. Here, we follow 
[TO] and use q = 5 with Kingslets. A discussion of 
the effect of varying q would be of interest but is 
out of the scope of the present work. 

The added complexity of running this algo- 
rithm, as compared to the standard Fourier- 
Galerkin method, comes from the forward and in- 
verse Fourier and wavelet transforms, and the iter- 
ations required to determine the threshold. Since 
the standard 4-th order Runge-Kutta scheme al- 
ready requires 12 Fourier transforms per timestep, 
the additional Fourier transforms represent an in- 
crease of computational cost of about 17% in total. 
The cost of each wavelet transform is proportional 
to Slog 2 (N) where S is the length of the wavelet 
filters, and for efficient implementations it is lower 
than the cost of a Fourier transform. Finally, the 
cost of the iterations is more difficult to evaluate 
since their number is not known a priori, but we 
observe in practice that it is low compared to the 
other costs. 




In Fig. 2a we show the solutions a few time steps 
after the appearance of the resonances, which do 
not occur for the CVS-filtered solution (shown in 
black). Figures 2b and 2c show that the evolu- 
tion is stable and we still have no trace of reso- 



Figure 3: Zoom of resonance at t = 0.037. Green 
(gray): Galerkin-truncated solutions. Black: CVS- 
filtered with Kingslets. 



nances, even for longer integration times when the 
Galerkin-truncated solution becomes perturbed, al- 
though after the formation of shocks the Gibbs phe- 
nomenon is intense (as discussed in [10] [9]). In 
Fig. [3] we show in detail how the resonances are 
completely filtered out by the CVS method. 

To demonstrate that the whole dynamics of the 
Burgers equation is preserved by CVS filtering, we 
plot in Fig. [4] the filtered profile along with the 
analytical solution as a reference, calculated using 
a Lagrangian map method [17]. One sees a very 
good agreement with only small discrepancies at 
the shocks due to the Gibbs phenomenon. 

Overall it appears that this implementation of 
the CVS filtering method achieves sufficient energy 
dissipation at shock locations to keep the numer- 
ical solution close to the desired entropy solution. 
It would be interesting to understand which ele- 
ment in the algorithm is essential for this benefi- 
cial dissipative effect, but unfortunately there are 
several competing influences which are difficult to 
disentangle. The filtering operation in itself (dis- 
carding the incoherent coefficients) is certainly an 
important source of dissipation, but it is difficult to 
quantify a priori since the Kingslets complex- valued 
wavelets are not an orthogonal basis, but merely 
a tight frame (see [7 ). Moreover, the alternating 
projections between the Fourier basis and a wavelet 
basis, which do not commute which each other, also 
introduce some dissipation. A first step in order to 
better understand the process by which this filter 
achieves dissipation is to move from a wavelet frame 
to an orthogonal wavelet basis, as we discuss in the 
next section. 

2.3 Elimination of resonances using 
real-valued orthogonal wavelets 

Although the Kingslet frame is well suited to sup- 
press resonances as we have recalled in the pre- 
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Figure 4: Solutions of the truncated inviscid Burgers equation at t = 0.037 (a), t = 0.048 (b), and 
t = 0.129 (c). Green (gray): CVS filtered with Kingslets. Black: analytical solution. 



vious section, it is appealing to be able to use a 
non-redundant real-valued orthogonal wavelet ba- 
sis. Due to its lack of translation invariance, this 
kind of basis does not perform well in the context 
of the algorithm described in the previous section 
[9] . Following previous work on CVS filtering of the 
2D and 3D Navier-Stokes equations [4j[T5j[Tl], we 
introduce the concept of a safety zone in wavelet 
space, that is, after computing the coherent coef- 
ficients as in the 4 th step of the CVS algorithm, 
we also keep the neighboring wavelet coefficients in 
space and in scale. The aim is to account for trans- 
lation of shocks to neighboring positions and gen- 
eration of finer scale structures from coarser ones. 
Hence, we have to add a step 4b. to the algorithm 
described in section [2?2l as follows: 



4b. Definition of the safety zone in wavelet space. 
We create an index set A containing pairs 
A = (j, i) indexing each coherent wavelet co- 
efficient in scale j and position z, kept in step 
4. We then define an expanded index set A* 
including the neighboring coefficients in posi- 
tion and scale, namely, for each pair (j, i), the 
pairs depicted in Fig. [5] [14 . Finally, all the 
coefficients not present in A* are set to zero. 

This additional step is able to generate a more sta- 
ble solution, but the fluctuation level is still high 
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Figure 5: Definition of the safety zone around the 
point (j,i). 



when compared to Kingslets. In order to smooth 
out these fluctuations we need a higher threshold in 
the CVS filter step of the algorithm, so we choose 
q = 8 in equation Q, changing accordingly the 
start-up value To. 

As examples we employ two different wavelet 
bases that are widely available in numerical analy- 
sis packages, the Daubechies 12 wavelet, which has 
compact support, and the Spline 6 wavelet, which 
has an exponential decay [2] . If we simply apply the 
CVS filtering procedure from section |^2] with these 
bases, the solution becomes very oscillating as soon 
as the resonances appear and we end up with poor 
results. But once the safety zone in wavelet space 
is implemented as described above, the dynamics is 
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Figure 6: Solutions of the truncated inviscid Burgers equation, CVS-filtered with a periodic Daubechies 
12 basis, at t = 0.067 (a), t = 0.143 (b), and t = 0.200 (c). Green (light gray): no safety zone. Red (dark 
gray): with safety zone. Black: analytical solution. 
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Figure 7: Solutions of the truncated inviscid Burgers equation, CVS-filtered with a periodic Spline 6 
basis, at t = 0.067 (a), t = 0.143 (b), and t = 0.200 (c). Green (light gray): no safety zone. Red (dark 
gray): with safety zone. Black: analytical solution. 
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Figure 8: Time evolution of energy. Energies of the 
CVS filtered solutions with safety zone collapse to 
the analytical energy evolution. 

properly preserved. In Figs. [6] (Daubechies 12) and 
(Spline 6) we see the significant improvement in 
the filtering capability of the code, comparing the 
cases with and without safety zone along with the 
analytical solution. 

The naturally oscillating character of real- valued 
wavelets and their lack of translation invariance still 
plays a role generating small perturbations (that 
get worse next to regions affected by the Gibbs phe- 
nomenon). But while the dynamics is lost when 
there is no safety zone, with huge oscillations cor- 
rupting the phase coherence of the shocks, after the 
introduction of the safety zone it is very well pre- 
served. Considering the time evolution of energy 
(Fig. [8]), it appears that in absence of a safety zone, 
not all the necessary energy is dissipated. This 
could be an explanation for the poor performance 
of the filtering scheme in that case. 

To give a quantitative aspect to the idea of "good 
filtering" we consider the global energy error esti- 
mate 



Jo v re f(x) 2 dx 



(5) 



where v re f is the reference analytical solution. This 
allows us not only to evaluate how close to the refer- 
ence the CVS-filtered solutions are, but also to com- 
pare the efficiencies of different wavelet bases. In 



Fig. 9a we plot the time evolution of e for all runs. 
For the unfiltered Galerkin-truncated solution, the 
error grows very fast as soon as resonances appear. 
The growth is slower for the CVS solutions without 
safety zone, but the solutions are still eventually 
destroyed. Due to their much smaller values, the 
error estimates for the Kingslets and for the real- 
valued wavelets with safety zone are barely seen in 
this plot. So in Fig. [9b] we change scales to find 
that they are of the order 10 -3 and stabilize once 
the influence of the resonances has been damped. 
We see that the errors of the real-valued orthog- 
onal wavelets stabilize very close to the Kingslets 
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Figure 10: Evolution of the fraction of incoherent 
wavelet coefficients. 



value. This makes their use attractive, a fact even 
more reinforced when we compare the level of com- 
pression along the time evolution (Fig. 10), i.e., the 



percentage of discarded coefficients each time step. 

Indeed, during a large fraction of the evolution, 
Kingslet-based CVS filtering keeps many more co- 
efficients than its counterparts based on orthogonal 
wavelets. The level of compression tends to stabi- 
lize at a slightly smaller value than the average of 
the other cases, but since the Kingslets frame has 
twice as many coefficients as real- valued orthogonal 
wavelet bases, this result shows the strong potential 
of the latter for the development of fully adaptive 
methods, provided a safety zone is implemented. 

3 2D Euler equation 

The emergence of resonances in Galerkin-truncated 
solutions of the 2D Euler equation was also shown 
in [12] . The fact that the CVS solutions, filtered 
with a 2D version of the Kingslets, are similar to 
the ones obtained from 2D Navier-Stokes with small 
viscosity [9] suggests that CVS might be suitable 
to filter the resonances in this case as well. There- 
fore, in the same spirit as in section |Z2} we apply 
the CVS method using Kingslets to the same initial 
condition used in the 2D example of [12]: 



2\k\ 



7/2 



-fc 2 /4+i0 k 



(6) 



where #k is a realization of a random variable uni- 
formly distributed in [0, 27r], k is the integer part of 
|k|, and Nk is the number of distinct vectors k such 
that fc < |k| < fc+1. The particular realization used 
in [12] as well as here can be retrieved online]^] The 
2D Galerkin-truncated Euler equations are solved 
using a fully dealiased pseudo-spectral method at 
resolution N 2 — 1024 2 with a low storage third 
order Runge-Kutta scheme for time discretization. 



A http : //www . kyoryu . scphys . kyoto-u .ac.jp/ 
°/ 7Etakeshi/populated 




Figure 9: (a) Time evolution of the global energy error e. (b) Rescaling of (a). 



The time step is adjusted dynamically to satisfy the 
CFL stability criterion. For more details on the nu- 
merical method, we refer the reader to [9]. 

In contrast to the Burgers case previously pre- 
sented, we do not have here an analytical solution 
to compare with, and make an error estimate, but 
a visual qualitative comparison will be sufficient to 
check if CVS filters out the resonances while pre- 
serving the dynamics. Resonances are well exhib- 
ited in plots of the Laplacian of vorticity, so, fol- 
lowing [12], we show contours of this quantity at 
t = 0.71. Figure [TT] shows the contours for the 
whole domain and we can easily see that CVS solu- 
tions do not show the resonances but keep the same 
general aspect. A more precise comparison can be 



made from Fig. 12 , where the contours of both cases 
at t — 0.71, zoomed-in around a region of intense 
resonance, are plotted together (left panel), as well 
as a cut as a function of x\ along a segment near 
x 2 = 3 (right panel). One sees very well how the 
resonances are suppressed and how the profiles are 
strikingly similar, indicating that the filter is able 
to maintain the physical aspects of the solutions. 
Finally, the dissipative character of the CVS filter 
is confirmed when considering the time evolution of 
the enstrophy Z = \ Jcj 2 , as shown in Fig. 
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4 Conclusion 

The continuous wavelet transform allowed to get 
further insight into the scale-space dynamics of res- 
onance phenomena in Galerkin truncated inviscid 
equations. We showed that oscillations appear in 
a non-local fashion as soon as a shock affects the 
cut-off scale, and that the resonant points and the 
shock act as sources of perturbations at the cut-off 
scale. We could also see that despite the fact that 
the resonances first appear at small scales, large 
scale structures develop at the resonant points and 
are stretched into smaller scales until they reach 
the cut-off and start acting as new sources of trun- 
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Figure 13: Time evolution of enstrophy for the 
Galerkin-truncated and CVS-filtered Euler equa- 
tions. 



cation waves. These new perturbations spread and 
reach the shocks, leading to energy equipartition. 

For the ID inviscid Burgers equation, the re- 
sults presented here confirm that the CVS filter- 
ing method we have previously proposed in [5] , us- 
ing a dual-tree complex wavelet frame (Kingslets), 
is well suited for eliminating all spurious oscilla- 
tions present in the Galerkin-truncated solution as 
reported in [12 . The resonances, which are not 
due to the dynamics of the original equation but 
rather to its discretization by a Galerkin method, 
are completely suppressed in this approach. Their 
'incoherent' character in relation to the system evo- 
lution is established. In order to better understand 
the dissipative process characteristic of CVS filter- 
ing, we have sought to replace Kingslets by stan- 
dard real- valued orthogonal wavelets. We have ob- 
tained satisfactory results under the condition that 
the coefficients which are adjacent to those whose 
moduli are above the threshold value are preserved. 
Such a safety zone is only necessary with orthogo- 
nal wavelets, to compensate for their lack of trans- 
lation invariance, as originally introduced for CVS 
filtering of the 2D and 3D Navier- Stokes equations 

mug eu. 

For the 2D Euler equation we have shown that 
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(a) (b) 

Figure 11: Contours of the Laplacian of vorticity at t = 0.71 (from —200 to 200, increments of 25). (a) 
Galerkin-truncated. (b) CVS. 




Figure 12: (a) Zoomed contours of the Laplacian of vorticity at t = 0.71. Gray: Galerkin-truncated. 
Black: CVS. (b) Laplacian of vorticity at t = 0.71 along a segment near to #2 =3. Gray: Galerkin- 
truncated. Black: CVS. 
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CVS filtering with Kingslets is also capable of fil- 
tering the resonances without perturbing the dy- 
namics. The filtered solutions match the unflltered 
ones but for the non-physical oscillations which are 
eliminated. The authors of [12] asserted that many 
features of the resonance phenomena were also ob- 
served in the 3D Galerkin-truncated Euler equa- 
tions, though these results have not been reported 
yet. It is an interesting perspective to test if in 
this case CVS filtering is still able to eliminate the 
resonances. 

A limitation of the approach presented here is 
that the solution is transformed back and forth at 
each timestep between the wavelet and the Fourier 
truncations, which do not commute with each 
other. These alternating projections are likely to 
introduce a weak dissipation in addition to the fil- 
tering operation per se. Therefore from the present 
results it cannot yet be determined whether the ob- 
served elimination of resonances could be achieved 
solely with wavelet filtering, or whether the inter- 
leaved truncations in Fourier space play a crucial 
role. This question could be answered by apply- 
ing the filtering method to the Wavelet- Galerkin 
truncation of the equations, instead of the Fourier- 
Galerkin truncation that was considered here, of- 
fering an appealing perspective for future work. 
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